Inequalities for Light Nuclei in the Wigner Symmetry Limit 
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Abstract 

Using effective field theory we derive inequalities for light nuclei in the Wigner symmetry limit. 
This is the limit where isospin and spin degrees of freedom can be interchanged. We prove that the 
energy of any three-nucleon state is bounded below by the average energy of the lowest two-nucleon 
and four-nucleon states. We show how this is modified by lowest-order terms breaking Wigner 
symmetry and prove general energy convexity results for SU{N). We also discuss the inclusion of 
Wigner-symmetric three and four-nucleon force terms. 
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Weinberg was the first to apply effective field theory to the nucleon-nucleon interaction 
Q]. Since then effective field theory methods for nuclear physics have been further developed 
and applied successfully to the two and three-nucleon system P S] P Recently though 
there has also been progress in applying effective field theory to many-body nuclear physics 

0030 

The traditional approach to many-body nuclear physics is based on two and three-body 
potential models. Although this approach has been very successful in describing the ground- 
state properties of light nuclei and neutron drops, the effective field theory approach offers 
several important advantages. First of all effective field theory provides a direct connec- 
tion with quantum chromodynamics (QCD). It explains the form of the interactions and 
how the strength of the interaction changes with the cutoff scale. Also the calculations 
can be systematically improved and one can estimate the errors due to contributions that 
have been neglected. Furthermore it has been discovered that the effective Lagrangian is 
often amenable to analysis so that non-perturbative results can be proven without even any 
numerical simulations. Effective field theory was used in to prove inequalities for the 
correlations of two-nucleon operators in symmetric nuclear matter. In upper bounds 
were proven for pressure in isospin-asymmetric nuclear matter and neutron matter in a 
magnetic field. 

The idea of proving rigorous inequalities from Euclidean functional integrals is not 
new. In fact the connection between effective field theory inequalities and nuclear 
attice simulations is much the same as the connection between QCD inequalities 



QHHIitIHQQQ 



and lattice QCD. In this letter we first confirm the work 
22| by showing that the physics of low-energy symmetric nuclear matter is close to the 
Wigner limit, where the isospin and spin degrees of freedom can be interchanged. In 
this limit SU{2) x SU{2) spin-isospin symmetry is elevated to an SU{4) symmetry. We 
prove that the energy of any three-nucleon state is bounded below by the average energy of 
the lowest two-nucleon and lowest four-nucleon states. We show how this is modified by 
lowest-order terms breaking Wigner symmetry and prove general energy convexity results 
for SU (N). We also discuss the inclusion of Wigner-symmetric three and four-nucleon force 
terms. Much of the physics we describe is universal and appears in other systems such as 
trapped Fermi gases near a Feshbach resonance 2^. Our general result for SU{N) will 
be interesting if and when one is able to trap four or more degenerate fermionic states. 



2 



Let represent the nucleon fields. We use r to represent Pauli matrices acting in 
isospin space and a to represent Pauli matrices acting in spin space. We assume exact 
isospin symmetry. In the non-relativistic limit and below the threshold for pion production, 
we can write the lowest order terms in the effective Lagrangian as 



^ = ^[^^0 + 2S7 - « - Z^)]^ - \CsNNNN 
- \Codd [NaN ■ NaN - NtN ■ NrN] . 



(1) 



We have neglected three-nucleon terms for now but will consider them later. We have 
written the Lagrangian so that the operator multiplying Codd Aips sign under the exchange 
of isospin and spin degrees of freedom. 

We now calculate Cs and Codd on a spatial lattice for various lattice spacings aiatuce- 
We sum all nucleon-nucleon scattering bubble diagrams on the lattice, locate the pole in 
the scatterin g am plitude, and compare with Liischer's formula for energy levels in a finite 
periodic box j25|j26|. In Table 1 we have computed these coefficients. 

Table 1: Contact potential coefficients 



«raLce(MeV) 


Cs (MeV~^) 


Codd (MeV-') 


20 


-3.40 X 10"^ 


-3.8 X 10-5 


40 


-1.20 X 10"^ 


-6. X 10-*^ 


60 


-7.70 X 10"^ 


-2.4 X lO^"^ 


80 


-5.60 X 10-5 


-1.3 X 10-*^ 


100 


-4.40 X 10-5 


-8. X 10-6 



As noted in j22|j2S|, deeply-bound Efimov states will begin to appear if we make aj'^^i^e too 
large. In reality this is not much of a restriction since we should not let a~[attice exceed the 
chiral symmetry breaking scale (or even ttt-tt for the pionless theory). 

We see that Codd is much smaller in magnitude than Cs- In the limit Codd ^ the 
SU{2) X SU (2) spin-isospin symmetry is elevated to an SU (4) symmetry. This symmetry 
was first studied by Wigner j2^0]j^, and arises naturally in the limit of large number of 
colors, Nc Although the ^Sq and ^5"! scattering lengths are quite different, the fact 

that both scattering lengths are large suggests we are close to the Wigner limit 

When Codd = the grand canonical partition function is given by 



DNDNexp{-SE) 



DNDNexp / cI^xCe 



(2) 



where 

Ce = -^[^4 - + - ii)\N - \CsNNNN. (3) 
Using Hubbard-Stratonovich transformations, we can rewrite Zq as 

oc y DNDNDfexp (^J d^^xC^j > (4) 



where 

4 = -N[d4 - 2S7 + « - Z^)]^ + CsfNN + \Csf. (5) 

Since C5 < 0, the / integration is convergent. 

Let M be the nucleon matrix. M has the block diagonal form, 

M = Mbiock ® Mbiock ® Mbiock ® Mbiock, (6) 
where we have one block for each of the four nucleon states and 



-^block — 



^4 - + « - ^')\ + ^sf. (7) 



We note that M is real valued and therefore det M > 0. 

Consider the two-nucleon operator A2{x) — [N\i[N\j{x), where i ^ j are indices for two 
different nucleon states. The two-point correlation function for A2 is 



Using our Euclidean functional integral representation, we have 

{A2{x)A\{Q)) = I DQ [M^,l^{x,Q)\' (9) 

where DQ is the positive normalized measure defined by 

Df dct M exp {\Cs [ d^x p) , , 

DO = — ' (10) 

/ Dj det M exp (|Cs / rf^^ /2) ' ^ ^ 

We note that since Mbiock is real valued, M^^^j^^ is also real valued. 

Next we consider the three-nucleon and four-nucleon operators A-i{x) — [-^]i[-/V]j[A'"]fc(x) 
and ^4(2;) = where i, j, A;, I are all distinct. We have 

A3(a;)4(0)) = \r>Q [MbIo^k(^,0)]^ (11) 

^4(x)4(0)) = {M^U{x,^t- (12) 



We note that 



-y [^4lock(^,0)]y j DQ [M^,Ui^,0)]\ (13) 
where the second hne is from the Cauchy-Schwarz inequahty. Therefore 

%{x)Al{0)) ^ < J(A,{x)Al{0))^^{A,{x)Ai{0))^^. (14) 

Let be the energy of the lowest state that couples to A2, and E^^ be the energy of 
the lowest state that couples to A4. Taking the limit x — > cxd in the temporal direction we 
conclude that any state with the quantum numbers of A3 must have energy less than the 
average of and Ea^, 

Eas > \ [Ea, + Ea,] . (15) 

Taking the limit x ^ 00 in any spatial direction we have that the inverse correlation length 
for A-i must be greater than the average of the inverse correlation lengths for A2 and A^, 

Ca] > I [CaI + CaI] . (16) 

Using arguments similar to those in we can show that the inequalities ()15j) and ()lfi|l 
hold for a general three nucleon operator 

u 

n 



^3(2;) = / d xid X2aijk{xi,X2)[N]i{x + Xi)[N]j{x + X2)[N]k{x + X3), (17) 



so long as fl is bounded. 

In the real world Codd is small but nonzero. We can measure the shift in the energy of a 
given state 1^4) using first-order perturbation theory, AEa = {A\ H' \ A), where 

H' = \Codd J d^x [NaN ■ NaN - NtN ■ NtN] . (18) 

Let us first consider two-nucleon states in an S-wave. In terms of SU{4) representa- 
tions, we decompose the tensor product of two fundamental 4 dimensional representations, 
4®4 = 10 ©6. The 6 dimensional representation is antisymmetric, and the spin and isospin 
representations must be 6 = (1, 0) © (0, 1). The spin triplet with isospin singlet corresponds 
with the deuteron, D, while the spin singlet with isospin triplet corresponds with the nearly 



bound ^5*0 states. A similar analysis for 5- wave three-nucleon states gives us one antisym- 
metric 4 representation with spin-isospin content 4 = (|,|)- This multiplet corresponds 
with the triton, T, and Helium-3. There is only one antisymmetric iS-wave four-nucleon 
state. It is a therefore spin singlet and isospin singlet and corresponds with Helium-4. 

Under a transformation that interchanges spin and isospin degrees of freedom, |^He) is 
mapped into itself, possibly with a minus sign. However H' is odd under this transformation 
and therefore, (^He| H' |^He) = 0. Under the interchange of spin and isospin, the spin-up 
^He state and spin-down T state are also mapped into themselves, possibly with minus 
signs. We conclude that (^He| H' |^He) = (T| H' \T) = 0. Under the interchange of spin 
and isospin, the D states and ^Sq states interchange with each other, again possibly with 
minus signs. Therefore {D\ H' \D) = — {^So\ H' \^Sq). We can now adjust for the first-order 
energy corrections due to H' , 



-EsHe, Et>\\_\ {Ed + EiSq) + i?4He] • 



(19) 



The physical binding energies are shown in Table 2 

Table 2: Binding energies for light nuclides 





~ MeV (nearly bound) 


D 


-2.224 MeV 


3He 


-7.718 MeV 


T 


-8.481 MeV 


^He 


-28.296 MeV 



Plugging these values into (fTIH) . we find that the inequality is satisfied, —7.7 MeV, —8.5 
MeV> —14.7 MeV. An analogous relation can be derived for the inverse correlation lengths. 



?3He' ST — 



(20) 



We stress that all of these inequalities hold in symmetric nuclear matter at any density 
or temperature where the effective theory description is valid. From here on we will con- 
sider only energy inequalities since the inverse correlation length inequalities are completely 
analogous. 

We now generalize the results ()15p and p6|) to the case where the Wigner symmetry is 
an SU{N) symmetry, for arbitrary A^. Let ngmaiunug, and n be any integers such that 



6 



< 2nsmaii <n < 2nug < N. Then 



E^> ^""^"^"^ E2n,- + ^""'^^ E2n „• (21) 



This inequahty is a statement of convexity of the energy as a function of nucleon number, 
with the additional requirements that the number of nucleons is less than or equal to 
and the two endpoints have an even number of nucleons. The proof of the inequality is 
straightforward. We can write 

^0 |^block(^,0)r = J DQ \M-,Ux,0)\'\M-l^{x,0)r\ (22) 

Applying the Holder inequality to the right-hand side, one finds the upper bound (j^lj) . 

Taking = 4 and setting 2nsmau = 2, n = 3, 2nbig = 4, we recover (fTH|) . If however we 
let 2nsmaii = 0, n = 3, 2nhig = 4, we get > |i?4. There are no first order H' corrections 
in this case, and we see that in the real world this inequality is satisfied, —7.7 MeV, —8.5 
MeV > -21.2 MeV. Setting 2nsmau = 0, n = 2, 2nbig = 4, we get E2 > ^E^. With first 
order H' corrections we have 

I {Ed + El So) > |-E4He, (23) 

and this is also satisfied, -1.1 MeV > -14.1 MeV. 

Up to this point we have ignored three and four-nucleon forces. It has been shown that the 

nn 

dominant three-nucleon force is Wigner-symmetric j27[[28[. We now show that introducing 
Wigner-symmetric three and four-nucleon forces do not spoil positivity of the Euclidean 
functional integral so long as the three-nucleon force is not too strong and the four-nucleon 
force is not too repulsive. We want to find a Hubbard- Stratonovich transformation that 
reproduces a contribution to the action of the form, 

JJexp [c2[NN{x)]^ + C3[NN{x)f + c^NNix)]^] . (24) 

X 

Let us just concentrate on what happens at a single point x, and in our notation we suppress 
writing the x explicitly. We note that NN raised to any power greater than 4 mush vanish. 
So we have 

exp L {NnY + C3 {NNf + C4 {NnY 

= ao + a^NN + f {NnY + f {NNf + ^ {NnY , (25) 
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Co = 1, ai = 0, a2 = 2c2, as = 6C3, 
We now to try find a real function g{f) sucli tliat 



04 = 12c2 + 24c4. 



(26) 



/oo 
# exp [/iViV + gif)] =ao + aiNN + f (iViV)' + f (iViV)' + ff (iViV)' . (27) 
-00 

We observe tliat a Hubbard-Stratonovicli transformation of this form maintains the positive 
functional integral measure. Expanding the left-hand side, we have 



ao 


ai 


a2 




1 





2C2 


ai 


a2 


as 







2C2 


6cs 


02 


as 


a4 




2C2 


6cs 


124 + 24c4 



a„,= / d/r exp [(?(/)], n = 0,1, 2, 3, 4. (28) 

Finding sufficient and necessary conditions for the existence of g{f) is known in the math- 
ematics literature as the truncated Hamburger moment problem. This problem has been 
solved jsSjjs^, and in our case g{f) exists if and only if the so-called block-Hankel matrix. 



(29) 



is positive semi-definite, with the added condition that if C2 = then C4 = 0. The determi- 
nant of this matrix is I6C2 — 8603 + 48C2C4. With an attractive two-nucleon force and small 
three and four-nucleon forces, the conditions are clearly satisfied. Whether or not these 
conditions are satisfied in the real world and at which lattice spacings is beyond the scope 
of this letter. But hopefully this will be numerically determined in the near future. 

To summarize, the physics of low-energy symmetric nuclear matter is close to the Wigner 
limit. We have proven that the energy of any three- nucleon state is bounded below by the 
average energy of the lowest two-nucleon and four-nucleon states. We have calculated the 
corrections due to the lowest-order terms breaking Wigner symmetry and shown that the 
inequalities are satisfied. We have proven general energy convexity results for SU{N) and 
shown that all of these inequalities are satisfied for = 4. We have also discussed the 
inclusion of Wigner-symmetric three and four-nucleon forces. 

The authors thank John Thomas for discussions and Paulo Bedaque for his hospitality 
at Lawrence Berkeley Laboratory where part of this research was completed. This work 
was supported by the U.S. Department of Energy under grants DE-FG-88ER40388 and 
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